Quantum-mechanical calculation of Stark widths of Ne VII n — 3, An = transitions 
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The Stark widths of the Ne VII 2s3s — 2s3p singlet and triplet lines are calculated in the impact 
approximation using quantum-mechanical Convergent Close- Coupling and Coulomb-Born-Exchange 
approximations. It is shown that the contribution from inelastic collisions to the line widths exceeds 
the elastic width contribution by about an order of magnitude. Comparison with the line widths 
CIh' measured in a hot dense plasma of a gas-liner pinch indicates a significant difference which may be 

naturally explained by non-thermal Doppler effects from persistent implosion velocities or turbulence 
developed during the pinch implosion. Contributions to the line width from different partial waves 
and types of interactions are discussed as well. 
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I. INTRODUCTION 



Spectral line shapes can provide very rich and valuable information on important plasma parameters, such as 
ion and atom temperature, electron density, electric field distributions, etc. The quantum-mechanical theory of 
collisional impact line broadening is well established and developed jy, however, the number of purely quantum 
[ — ■ calculations, especially for highly charged ions, is rather limited. Most theoretical efforts were directed toward 
I/"") ' elaboration of semiclassical or semiempirical methods which showed good accuracy for neutrals and low-charge ions. 
It is only recently, when a number of sophisticated atomic collisional codes have become available, that high-quality 
quantum-mechanical results could be applied to line shape calculations for highly charged ions. From the experimental 
point of view, test measurements of line profiles are impeded by the required independent determination of plasma 
temperature and density. The experimental situation is even more peculiar in that the line widths of high-Z ions were 
measured almost exclusively by the Bochum group (see and references therein), and therefore lack an independent 
Y\ • confirmation. 

The recent results on the Stark broadening of spectral lines from multiply charged ions revealed a significant discrep- 
ancy between the independent quantum-mechanical calculations and, on the other hand, experimental measurements 
and semiclassical results. For the Bill measurements the Stark line widths for the simplest 2s — 2p transition 
O r differ by as much as a factor of 2, the two quantum results |||7| being in agreement to within 10%. A possible 
explanation for this discrepancy in terms of a developed turbulence and different treatments of small partial waves in 
electron-ion scattering was proposed in Ref. Q; however, more comparisons and detailed investigation of important 
contributions to the line width are of primary importance. 

Measurements of line profiles for the 2s3s — 2s3p transitions 1 Sq — 1 Pi and 3 Si — 3 P2 of Ne VII emitted from a hot 
dense plasma of a gas-liner pinch were reported recently ^,pf- The experimental line widths for singlet (A s = 3643.6 
A) and (the strongest) triplet (A T = 1982.0 A) lines are AA^ = 1.70 ±0.26 A and AA T = 0.45 ±0.07 A, respectively. 
The electron density and temperature were measured independently by laser Thomson scattering and turned out to 
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be in the ranges N e = (3 — 3.5) • 10 18 cm -3 and T e = (19 — 20.5) eV. The measured line widths agree well with most 
semiclassical [l]]^,|§] or semiempirical [fiof calculations but only marginally with other semiempirical results [fll| . 

Here we present the results of fully quantum-mechanical calculations of the Stark line widths for the n = 3 — > 3 
transitions in Be-like Neon. The plan of this paper is as follows. In Sec. II the calculational method is described. The 
features of atomic structure as well as inelastic and elastic contributions to the line widths are discussed in detail. In 
Sec. Ill a comparison with available experimental and theoretical results is made and the sources of discrepancy are 
investigated. Finally, Sec. IV contains conclusions and recommendations. 



II. METHOD 




A. General theory 

The calculational method applied here is basically the same as the one used in Ref. ||. We start from the 
fundamental formula for the full collisional width at half-maximum (FWHM) for an isolated line corresponding to a 
transition u — > I [|l2]| : 

oc / 

w = N e I vF (v) ( * uu , (v) + £ a u , (v) + / \f u (0, v) - ft (0, v)\ 2 dQ \ dr. ( I ) 

i'^i 

with N e being the electron density, v the velocity of the scattered electron, and F (v) the Maxwellian electron velocity 
distribution. The electron impact cross sections a uu > {a u > ) represent contributions from transitions connecting the 
upper (lower) level with other perturbing levels indicated by primes. In Eq. (Q), the f u (0,v) and fi (9,v) are elastic 
scattering amplitudes for the target ion in the upper and lower states, respectively, and the integral is performed 
over the scattering angle 8, with dH, being the element of solid angle. Equation (|l|) relates a line width in the 
impact approximation with atomic cross sections, facilitating the use of well-developed techniques of atomic scattering 
calculations for line broadening studies. It can also be rewritten in terms of the elastic S-matrix elements assuming 
LS coupling (see, e.g., [^3[): 

[ LlLjS T ll' 



x [5,/, - S v (L u Sl'\L T u S T ; L u SI\lZs t ) S* l (USl 1 \hj S T ; L x Sl\h{ ' S T )\ } . (2) 

Here L and S are the atomic orbital angular momentum and spin, I and V are the electron orbital angular momentum 
before and after collision, superscript T denotes the quantum numbers of the total electron+ion system, and the 

^ g ^ 1 are 6-j symbols. The advantage of Eq. ([!]) is that it gives more clear insight as to the importance of 

inelastic and elastic contributions to the line width; therefore, we will mainly be referring to Eq. ([!]) in what follows. 

In the present work, the inelastic cross sections appearing in Eq. (0) were calculated with two independent methods, 
i.e., the Convergent Close- Coupling (CCC) and Coulomb-Born-Exchange (CBE) approximations. The basic idea of 
the CCC method lies in the close-coupling expansion with a large number of square-integrable states. A set of 
coupled Lippmann-Schwinger equations for the transition matrix is solved in momentum space, and the convergence 
of the results may be checked easily by increasing the number of the basis functions. The details of the CCC method 
can be found in a number of recent reviews [ |15|flq| , where a very good agreement with various experimental collisional 
data is shown. For the calculations with the more traditional CBE approximation we made use of the code ATOM 
described in Ref. ]17j . In addition to the Coulomb attraction between ion and electron and exchange, ATOM 
accounts for normalization (unitarization) effects and uses experimental level energies when calculating the atomic 
wave functions and collisional cross sections. It is well known that the Coulomb-Born approximation corresponds to 
perturbation theory with 1/Z as expansion parameter, where Z is the spectroscopic charge; therefore, one can expect 
better accuracy for the CBE method applied to highly charged ions. Although the CCC method generally provides 
a superior accuracy, the use of ATOM greatly reduces the computational efforts. Comparison of CCC and CBE 
cross sections for highly charged H- and Li-like ions (Z < 12) demonstrated an excellent agreement between these two 
methods and available experimental data [H. 
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B. Atomic Structure 



Before proceeding to the details of collisional calculations, it is worth mentioning some features of the Ne VII atomic 
structure. First, the current version of the CCC code utilizes the Hartree-Fock (HF) frozen-core approximation for 
atomic wave functions. To study the validity of this approach, we have made a comparison of wave functions calculated 
with the full HF and HF frozen-core methods using the Cowan code jnj . The agreement between the two sets of wave 
functions proved to be very good, thereby justifying the use of the frozen-core approximation. These calculations were 
also used to determine the root mean squared radii of the lower 2s3s and upper 2s3p states, which were found to be 
1.81ao and 1.79ao respectively, where ao = 0.529 • 10~ 8 cm is the Bohr radius. (Recall that the CBE ATOM code 
constructs atomic wave functions by solving the Schrodingcr equation with experimental energies rather than solving 
ab initio HF equations.) Another measure indicating the level of accuracy are the oscillator strengths, which for the 
CCC calculations are found to agree within one percent with the Opacity Project results (20). Since in some cases the 
ATOM oscillator strengths Jatom deviate from the high-accuracy results f acc by as much as 15%, the CBE dipole- 
allowed excitation cross sections were rescaled by the ratio face/ /atom to improve these results. Finally, to check the 
applicability of LS coupling, we carried out large scale atomic structure calculations for Ne VII with Cowan's code 
taking into account both intermediate coupling and configuration interaction. The results obtained show that the 
levels of interest of Ne VII correspond to practically pure LS coupling, although configuration interaction is important 
for the 2s3p 1 P and 2p3s 1 P states which mix to a level of 10%. Nonetheless, this mixture is unlikely to be important, 
since in the sum of inelastic cross sections in Eq. (fj]) this effect is essentially smoothed out. 



C. Inelastic collisions 

The inelastic cross sections appearing in Eq. (|l|) include all possible electron-induced transitions originating from 
the lower or upper states of a transition. It is normally safe to neglect the ionization and recombination processes, 
taking into account only electron impact excitation and deexcitation. For the line widths discussed here, even the 
An ^ excitations may (but will not) be ignored, since their rates are at least two orders of magnitude smaller 
than those for the An = transitions. Table I presents the CBE rate coefficients for electron impact excitation and 
de-excitation processes connecting the upper and lower levels of transitions with other perturbing 2131' levels. The 
calculation was carried out for an electron temperature T e = 20 eV, which corresponds to the experimental conditions 
of Ref. [0, and only one-electron transitions are considered here since two-electron transitions were found to have 
much smaller cross sections. One can see that the largest rate coefficients correspond to dipole-allowed transitions, 
while dipole-forbidden and spin-forbidden channels contribute only a few percent to the inelastic part of the line 
width. It should also be noted that since the reaction thresholds are smaller than 6-7 eV, the rate coefficients are 
rather insensitive to small (a few eV) variations in the electron temperature around the experimental value of 20 eV. 

Both calculational methods give close (within 10%) results for the most important dipole-allowed cross sections. 
(Note that the excitation of the inner 21 electron is also significant for the line width, contributing as much as 12% 
and 8% for singlet and triplet lines, respectively.) An example of the agreement between the CCC and CBE results 
is demonstrated on Fig. 1, where 2s3s 1 S - 2s3p 1 P and 2s3p 3 P - 2s3d 3 D excitation cross sections are shown. 
Unfortunately, there are no other available theoretical nor experimental data for the 3 — 3 transitions in Ne VII, so 
in order to test the accuracy of our calculations it seems to be reasonable to make a comparison with the existing 
2 — 3 data for this ion. Probably, the most accurate theoretical results were produced recently by Ramsbottom et 
al. pl[ , who calculated electron impact excitation rates for many 2 — 3 transitions using the multichannel R-matrix 
method. The comparison shows very good agreement between our data and those of Ref. f2l|| . For instance, the 
CBE excitation rate coefficients (T e = 10 6 K ss 86 eV) for the outer electron transition 2s2p 3 P - 2s3d 3 D and inner 
electron transition 2s2p 3 P - 2p3s 3 P are 6.8 x 10~ 10 cm 3 s~ 1 and 9.1 x 10~ n cm 3 s _1 , respectively, which agree well 
with the R-matrix values of 6.2 x 10~ 10 and 9.6 x 10~ n cm 3 s~ 1 . There also exist 0-pinch experimental results |22] ] 
for excitation rates from the ground and metastable states to some of the n = 3 states at an electron temperature 
of 260 eV; these are 2 to 3 times smaller than CCC/CBE rates, but large experimental errors up to 200-300% limit 
their usefulness. 

To summarize, for the experimental conditions of Ref. Jpj, the CBE inelastic contribution (with account of the 
An = 1 transitions) to the line widths obtained from Eq. (yTis wf n ~ 0.806 A for the singlet line and wf n « 0.197 A 
for the triplet line. 
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D. Elastic Collisions 



According to Eq. (|l|), the non-Coulomb elastic amplitudes of scattering from the upper and lower states at the same 
electron impact energy should be subtracted and averaged over the Maxwellian electron energy distribution. These 
amplitudes were calculated for a large range of electron energies only with the CCC code, since the existing version of 
the CBE code ATOM produces only inelastic cross sections. The 2s3s 3 S and 2s3p 3 P elastic cross sections a e i(E) 
along with the coherent difference term a(E) = J \f s (9, v) — f p (0, v)\ d£l are shown on Fig. 2, the singlet cross 
sections and difference term having a similar behavior. These results unambiguously reveal the same peculiarities as 
were noticed for the B III 2s — 2p elastic term ||, i.e., a faster than X/E energy dependence and strong cancellation 
in &(E). For example, at electron impact energies E > 30 eV, the coherent difference a is more than an order of 
magnitude smaller than any of the o~ e i. Since at large energies the elastic cross section is mainly determined by the 
size of a system, such a cancellation may be due to almost equal mean squared radii of the 2s3s and 2s3p states, as 
was already mentioned above. The general behavior of the elastic difference term deserves a special investigation and 
will be reported elsewhere. The Maxwell-averaged elastic contribution to the line width is u>f ; « 0.067 A and ~ 
0.023 A for singlet and triplet, respectively. This shows that in this case the elastic contribution to the line widths is 
about an order of magnitude smaller than the inelastic one, which is not surprising for such high temperatures. 



E. Final results 



To summarize, the total line widths (FWHM) for the 2s3s - 2s3p, 1 S - 1 P and 3 S - 3 P, transitions obtained from 
Eq. (P are wf » 0.873 A and wj « 0.220 A. The same widths were also calculated with Eq. (J|) using the CCC 
elastic T-matrix elements and the relation between T- matrix and S-matrix T = S — I (I is the unit matrix). The 
corresponding singlet and triplet widths are tuf ~ 1-05 A and w\ ~ 0.230 A. The difference between the results 
obtained with Eqs. ([!]) and (|^) can probably be attributed to the resonances in the CCC T-matrix, which were not 
included into the CBE inelastic calculations. A conservative estimate of the accuracy of these results, based on the 
CCC-CBE agreement and the accuracy of the CCC calculations along the [Be] sequence, is 15 %. Thus, the final 
Stark line widths are 

w s = (1.0 ± 0.15) A, w T = (0.23 ± 0.03) A. (3) 



III. DISCUSSION 

The line widths calculated here differ noticeably from the measured values of Ref. ^ and most of the theoretical 
data. The ratios of experimental to different theoretical Stark widths for the Ne VII lines are presented in Table II. The 
methods cited there cover various modifications of the semiclassical and semiempirical |I^,[Il] approximations. 
The semiclassical methods, including the latest nonperturbative calculations yield values which are generally in 
agreement with the experimental data. The semiempirical results of Dimitrijevic and Konjevic [Ol are rather close 
to our values, and this is quite similar to what had already been noticed for the B III calculations H. 

The major diagnostics challenge in the gas-liner pinch experiment § may be the determination of the main plasma 
parameters, i.e., the electron temperature and density, in a region where the multiply charged ions of neon are 
situated. In the experiment both T e and N e were determined from the Thomson scattering only globally which may 
not be characteristic of the plasma conditions near the locally injected neon. As a matter of fact, there exist some 
experimental indications that density and temperature do vary in the vicinity of the doping gas [|23fl . However, the 
experimental value of T e is supported by the fact that electron temperatures T e = 19 — 20 eV are well within the 
range of the maximal abundance temperatures for NeVII at an electron density N e = (3 4) • 10 18 cm~ 3 . Our 
calculations with the collisional-radiative code NOMAD |24|] show that for equilibrium conditions the Ne VII ions 
account for about 30% of the total amount of neon. Another line broadening mechanism affecting the observed widths 
may be unresolved Doppler line splitting associated with the radial implosion velocities in the gas-liner pinch |p3| . 
The contribution from an ion (proton) collisional broadening may be estimated using Eq. (517) from Ref. 0, and it 
is easy to show that ion broadening is negligibly small comparing to electron impact broadening. 

Since the experimental conditions in the NeVII measurements were basically the same as for the B III experiment, 
the general conclusions |(| regarding a possible effect of a developed turbulence on the line widths should remain 
essentially the same. It was mentioned in Ref. H that the measured value of Stark width for the triplet transition 
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AA w 0.45 A constitutes about 70% of the total measured line widthQ which therefore is AA cxp ~ 0.64 A. This full 
width includes Stark, Doppler and instrumental broadening, the latter being decomposed into Gaussian (0.07 A) and 
Lorentzian (0.05 A) parts B. For an ion temperature of T. t — T e — 20 eV, the pure Doppler width is approximately 
AAc w 0.15 A. As noted in H, the Reynolds numbers for the Bochum gas-liner pinch experiment are of the order of 
10 4 , which is sufficient for a developed turbulence to exist. Such a turbulence leads to an extra chaotic motion of Ne 
ions with a characteristic velocity of the order of the proton thermal velocity v p . Hence, the full thermal+turbulent 
Doppler width becomes a factor \/20 + 1 » 4.6 larger (here 20 is the ratio of masses Mn £ /Mh) and is now AAd « 0.70 
A. Using Eq. (6) of Ref. ^5|, for the FWHM of a combined Voigt profile including Stark, thermal+turbulent Doppler 
and instrumental contributions, we get a value AA rs 0.85 A which is 30% higher than AA CX p- The main uncertainty 
in this calculation obviously comes from the turbulent contribution, which is rather sensitive to the value of the 
characteristic velocity. It is straightforward to show that reducing this velocity by one third only, i.e., multiplying the 
pure Doppler width by 3.1 instead of 4.6, one can exactly reproduce the experimental value of AA CX p- Thus, according 
to the hypothesis proposed in Ref. Q , reasonable values of characteristic turbulent velocities may naturally explain 
the observed difference in line widths. 

Regarding the discrepancy between the quantum and other theoretical calculations, the reader may wonder as to 
the source of such a difference. The crucial point is that unlike the quantum-mechanical methods, the semiclassical 
approaches have a natural limit of applicability arising from the Heisenberg uncertainty principle (see, e.g., J2(|). The 
criterion of applicability of the semiclassical calculations may be formulated p7fl as a requirement for the distance 
of the closest approach r m j n , rather than the impact parameter p, to be larger or at least of the same order than 
the corresponding de Broglie wavelength, A m i n = 2Trh/mv max . Using the angular momentum conservation it is 
straightforward to show that this is equivalent to the inequality: 

2tt < L. (4) 

Another limitation on impact parameters was introduced in order to avoid violations of unitarity jlj , but still assuming 
the long-range dipole interaction to remain valid. Again reformulated in terms of the distance of the closest approach 
the corresponding condition may be written as 

— < 1, (5) 



where r n is the excited state atomic radius. If this inequality is violated, both semiclassical and long-range interaction 
approximation are questionable. Using the Coulomb parameter r\ = [Z — l)e 2 /hv, Eq. (Q) may also be expressed in 



terms of the total angular momentum L as 27 



[Z-D r„ [ 1 + iL/rjf 1 " 
a L 2 



<1 (6) 



As was noted above, the mean root squared radii of the 2s3s and 2s3p states are about 1.8ao. For T e = 20 eV, 
the Coulomb parameter is ij ~ 7 and therefore the ratio r n /r min takes values of 1.45, 0.96, 0.43 and 0.22 for L = 
4, 5, 8 and 12 respectively. It follows then that for the given electron temperature, criteria (|J) and (JsJ) are similarly 
restrictive for the semiclassical approximation. 

Unlike to the semiclassical method, in fully quantum-mechanical calculations the determination of the range of 
significant L-values is naturally accomplished by the partial wave expansion. In Fig. 3 the contribution of different 
total electron+ion angular momenta Lt to the CCC cross sections is shown for an incident electron impact energy 
of 20 eV for a number of transitions^. Naturally, the elastic cross sections are governed by the smallest values of 
Lt, which are concealed in the strong collision term of semiclassical calculations. The most important inelastic cross 
sections having the smallest thresholds reach 50% of their values only for L = 9 for which the l.h.s. ratio of Eq. 
(|J) is about 0.35. This number is probably already sufficiently small to justify the use of the long-range interaction 
approximation for L > 9; however, the restrictions following from Eq. (Q) are less obvious to have been overcome. 

Another discrepancy may come from other than dipole interactions. Although the monopole interaction was not 
explicitly included into nonperturbative semiclassical calculations the quadrupole transition 2s3s — 2s3g? was shown 
to account for about 15% of the line width. This value is in contradiction to the present results. As one can see from 



1 There is no information in Ref. B on the full line width of the singlet transition, so we will not discuss it in what follows. 



2 The CBE partial wave composition is practically the same. 
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Table I, the 2s3s — 2s3d quadrupole channel contributes only approximately 3% to the quantum-mechanical inelastic 
line width. If we take into account only those transitions that were considered in Ref . |9[ , then this number increases 
to 3.5%, still a factor of 4 smaller than the nonperturbative semiclassical result. 

These considerations clearly show that the accuracy of the semiclassical calculations may not be as high as it is 
often thought to be, and new calculations, both semiclassical and quantum-mechanical, are needed to better establish 
the limits of applicability for the non-quantum methods. 

IV. CONCLUSION 

A fully quantum-mechanical calculation of the Stark line widths for the singlet and triplet 2s3s — 2s3p lines of 
NeVII was carried out in the impact approximation with the use of accurate atomic data. Although the results 
obtained disagree with experimental and most theoretical results, a natural explanation for this disagreement can 
be suggested. On one hand, the measurements are not free from difficulties related to possible extra contributions 
from turbulence and unresolved Doppler shifts. This suggests an independent measurement of Stark widths of highly 
charged ions. On the other hand, the semiclassical calculations, not obviously producing accurate results for other 
than dipole interactions, may have problems when being applied to the small impact parameter region. In our opinion, 
the next important step in the development of Stark broadening theory would be a very detailed comparison between 
quantum and semiclassical results. 
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Figure Captions 

Fig. 1. Electron impact excitation cross sections for the transitions 2s3s 1 S* - 2s3p 1 P and 2s3p 3 P - 2s3d 3 D in 
Ne VII. CBE - dashed lines, CCC - solid circles. 

Fig. 2. Non-Coulomb elastic cross sections of Ne VII ions in 2s3s 3 S (solid line) and 2s3p 3 P (dot-dashed line) 
states, and the coherent difference term a(E) (diamonds). 

Fig. 3. Contribution of different total electron+ion angular momenta Lt to various elastic and inelastic cross 
sections. 



TABLE I. The CBE electron impact excitation and de-excitation rate coefficients for Ne VII in units of cm 3 s 1 for T c = 20 
eV. 





2s3s 11 S 


2s3s 'S 


2s3p X P 


2s3p d P 


2s3s 3 S 




7.46(-10) 


3.29(-10) 


6.10(-08) 


2s3s 1 S 


2.21(-10) 




5.17(-08) 


1.13(-10) 
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1.31(-07) 




5.30(-10) 


2s3p 3 P 


1.35(-07) 


8.44(-10) 
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2s3d 3 D 


9.96(-09) 


1.50(-09) 


1.48(-09) 


7.75(-08) 


2s3d 1 D 


4.10(-10) 
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6.45(-08) 
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2p3p X P 






1.20(-08) 


4.34(-ll) 


2p3p 3 D 






2.04(-10) 


9.69(-09) 


2p3p 3 S 






3.89(-ll) 


2.18(-09) 


2p3p 3 P 






1.12(-10) 


3.86(-09) 


2p3p 1 D 






7.04(-10) 


5.92(-ll) 


2p3p 1 S 






6.63(-10) 
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TABLE II. Ratio of the experimental Stark widths of the 2s3s — 2s3p lines in Ne VII to different theoretical widths. 



Line 




T e (eV) 


N c (cm" J ) 


w cxp /w theor 


Ne VII 'So 


- 1 Pi 


19 


3.5-10 18 


1.28" 


1.15" 


1.57 c 


0.88(0.77) d 


1.70 e 


Ne VII 3 Si 
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